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ABSTRACT 

We present a BIMA image of the supernova remnant G21.5— 0.9 at 94 GHz with angular resolution 
8.6" x 4.5". On scales larger than our synthesized beam, our results do not indicate any radial or local 
variations in the acceleration or synchrotron loss processes for the relativistic particles emitting at radio 
wavelengths. However, the FWHM size of the radio remnant is significantly greater than that of the 
X-ray remnant. Either any low- frequency spectral break is distributed across the remnant, or more 
probably the break frequency is in fact higher than ~ 100 GHz. 

Subject headings: ISM: individual (G21.5— 0.9) — radio continuum: ISM — stars: neutron — stars: 
winds — supernova remnants 



1. INTRODUCTION 

Pulsar-powered nebulae have a characteristic radio spec- 
trum; powering by the embedded neutron star, whether 
observed or not, creates a relatively flat synchrotron spec- 
, trum with a spectral index a (S v oc u a ) between and 
—0.4 at low frequencies. Above this regime one or more 
breaks occur in the spectrum as synchrotron losses, re- 
duced energy from the pulsar with time, and expansion of 
the nebula play a role. At high frequencies, synchrotron 
losses will dominate the spectrum and a break will oc- 
cur such that the higher frequency spectral index should 
be steeper by 0.5 (Kardashev 1962; Reynolds & Chanan 
1984). When the expansion, pulsar decay, and synchrotron 
losses are all important, an intermediate region with an 
intermediate slope will be present in the spectrum. The 
break points will generally diverge with time (Woltjer et 
al. 1997). In addition, non-thermal X-rays may come from 
the same or different parent populations of relativistic elec- 
trons. 

To determine the parameters of the Crab-like SNRs, it 
is necessary to measure the frequencies at which the spec- 
tral breaks occur and the actual spectral indices in the 
different spectral components. Also, because the expan- 
sion rate and the magnetic fields may vary with position 
within a remnant, it is important to image in the region of 
the break frequency to look for possible variations across 
the source. The lower break frequency generally falls some- 
where between 10 9 and 10 13 Hz (Salter et al. 1989b; Green 
& Scheuer 1992). The Crab itself has a break frequency 
of about 10 13 Hz, well above the mm band (Strom & 
Greidanus 1992). However, measuring spatial variations 
in the mm-wave region of the spectrum is a difficult task 
which has not been done for any Crab-like SNR, since 
the synchrotron emission is generally weak and, while the 
sources are often extended, high resolution observations 
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are needed to resolve structure across them. 

G21.5— 0.9 is one candidate which can be studied in 
detail. Integrated flux density measurements indicate a 
break frequency near 50 GHz and a total flux density of 
4 Jy at 84.2 GHz (Salter et al. 1989b). The object has been 
imaged with 7" x 4" resolution at both 4.9 GHz (Becker 
& Szymkowiak 1981) and 22.3 GHz (Fiirst et al. 1988). 
As expected, these two low-frequency images appear the 
same with a filled center and some structure. In addition, 
the center has very low polarization but is surrounded by 
a ring of highly polarized emission with a radial magnetic 
field. This may well indicate different conditions in the 
interior and around the outside of this SNR. 

Recent X-ray observations have painted a very different 
picture at higher energies (Slane et al. 2000; Warwick et 
al. 2001). Here the emission is seen over a greater area, 
and has several distinct components: an unpulsed, com- 
pact X-ray source, an extended core corresponding to the 
radio nebula, and a faint outer shell of radius about 150 
arcsec which is not observed in the radio. 

In this letter we present the first interferometric image 
above the previously suggested break frequency using the 
BIMA 1 array at 94 GHz. We compare this new image 
with images at 5 GHz, 22 GHz and 0.1-10 keV to look for 
structures corresponding to the spectral break seen in the 
integrated flux density. 

2. OBSERVATIONS 

Observations were obtained with the B, C, and D con- 
figurations of the BIMA array (Welch et al. 1996) between 
1998 October and 1999 March. To achieve suitable uv cov- 
erage, we used the technique of multi-frequency synthesis 
over both sidebands of the local oscillator. We measured 
two orthogonal polarizations to ensure accurate total in- 
tensities. In order to image fully the extended structure, 
and to assist with the recovery of low spatial frequency 
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information, we imaged the source with a 7-point hexag- 
onal mosaic. The primary beam of the 6.1-m antennas is 
Gaussian with FWHM 2.13' and 1.74' respectively at 90, 
and 110 GHz. 

The quasar 1743—038 was observed as a phase calibra- 
tor at 25 to 30 minute intervals, after 3 complete cycles 
of 7 pointings. The antenna gains were determined from 
observations of planets at short antenna separations; the 
uncertainty in the flux density scale is 10%. All data were 
obtained in an 800 MHz bandwidth using a digital correla- 
tor with 256 frequency channels. Single sideband system 
temperatures during the observations were between 170 
and 450 K (scaled to outside the atmosphere). 

Imaging was done with the MIRIAD software pack- 
age (Sault, Teuben, & Wright 1995) using a maximum 
entropy mosaicing algorithm (Sault, Staveley-Smith, & 
Brouw 1996), in which the multiple pointing centers were 
combined into a mosaiced image corrected for the pri- 
mary beam. To assess imaging reliability for a complicated 
source we must consider image fidelity (i.e. how well the 
image represents the real source brightness distribution), 
rather than the more familiar measures of thermal noise 
level or image dynamic range. To do this we re-calibrated 
and re-imaged G21.5— 0.9 many times by choosing subsets 
of the data, and using a variety of uv weighting schemes 
(natural, uniform and robust). The RMS difference be- 
tween these images was consistent with the thermal noise 
at 94 GHz (3 mJy beam -1 ). The peak differences be- 
tween these images was 12 mJy beam -1 , so differences 
between images corresponding to as much as 1.5 contours 
on figure 1 may not be significant. We use these empirical 
estimates of the RMS in evaluating the results. 

3. RESULTS 

Figure 1 shows the maximum entropy image with 
mean frequency 93.7 GHz. The convolving half-power 
beamwidth was 8.6" x 4.5". This image includes measure- 
ments on spatial scales up to 2 arcmin. The 94 GHz image 
is similar to the 5 GHz (Becker & Szymkowiak 1981) and 
22 GHz (Furst et al. 1988) images with comparable angu- 
lar resolution, but does not reproduce the detailed small- 
scale structure seen at 22 GHz, even after convolution to 
the resolution (8") used for detailed studies by Furst et al. 
(1988). If real, these features should be above the noise 
level in our image. It may be that the 'CLEAN' algo- 
rithm used by them produced a mottled appearance in the 
smooth emission (Cornwell, Braun, & Briggs 1999). The 
three interferometric images also differ somewhat from the 
single dish 90 GHz image (Salter et al. 1989a): the sec- 
ondary peak to the north east in that image is not seen 
with the interferometers. Given the difficulties inherent in 
single-dish imaging, we do not consider this discrepancy 
significant. 

The total flux density recovered in the final image is 
3.5 Jy measured within a radius of 45". The total flux 
densities measured on the individual images made from 
the 90 and 110 GHz data separately (3.4 Jy and 3.6 Jy) 
are consistent within the calibration errors of about 10%. 
Omitting the most compact array configuration reduces 
the measured total flux density, to 3.0 Jy, as expected. 
We conclude that the final image recovers all the emis- 

2 The X-ray data used were from HRC image #1406 (1999 October 25) 



sion corresponding to the known radio extent, with a total 
flux density 3.5 ± 0.4 Jy at 94 GHz, in agreement with 
the single-dish result at 91 GHz of 3.8 ± 0.4 Jy (Salter et 
al. 1989a). However, we note that the total flux density 
has been estimated in the deconvolution and has not been 
directly measured by the interferometer. 

3.1. Radial profile analysis 

In order to evaluate the spectral index distribution accu- 
rately, we must match the data as closely as possible, not 
only in angular resolution but also in the spatial sampling 
of the images. One possibility is to compare completely 
sampled images at each frequency. We have compared 
the normalized surface brightness in elliptical annuli of 
the 94 GHz, 22 GHz (Furst et al. 1988) and 0.1-10 keV 
images. 2 The 5 GHz maps (Becker & Szymkowiak 1981) 
were not available in digital format. The best-fit center for 
the analysis was determined empirically to be 18 h 33 m 33!4, 
-10°34'04" (J2000). The ellipses had position angle 45° 
and axial ratio 0.6. To make the comparison we subtracted 
base levels of 5000 K, —2 mJy beam -1 and 2 counts re- 
spectively from the 22 GHz, 94 GHz and X-ray images, 
and then normalized the data at the center, after remov- 
ing the contribution of the compact X-ray source. The 
comparison, using annuli of width 5 arcsec, is plotted in 
figure 2. 

The errors in our radial profile at 94 GHz were estimated 
empirically from the RMS variation between the various 
images discussed in §2. Using 1.5 arcsecond bins, the RMS 
variation between the derived radial profiles was generally 
less than 1.5% of full scale. The error bars on figure 2 
indicate this uncertainty. 

The two radio datasets agree within our estimate of the 
error at all but two points on the plot. Taking into account 
the sensitivities of the images and the 'peak' errors in the 
imaging process (§2), we do not find these discrepancies 
to be significant. Thus we can find no significant radial 
differences between the radio images at 22 and 94 GHz. A 
similar analysis was made from the contour maps of the 5 
and 94 GHz data using graphical techniques. These were 
also found not to differ significantly. 

By contrast, the X-ray shape of the plerion remnant dif- 
fers significantly from the radio shape. At a radius of 30", 
the radio surface brightness is two-thirds of the central 
value, whereas the X-ray surface brightness is only 15% of 
its peak. The major axis FWHM sizes of the radio and 
X-ray remnants are 68" and 35" respectively. 

3.2. Search for a radio shell 

G21.5— 0.9 has recently been shown to have faint X- 
ray emission extending well beyond the radio (Slane et al. 
2000; Warwick et al. 2001). This may be the signature 
of the SNR shell; in this case a radio counterpart would 
be expected. Slane et al. report that they have been un- 
able to detect this shell at 1 GHz to a surface brightness 
limit of 4 x 10~ 21 W m~ 2 Hz -1 sr _1 (30 mJy arcmin -1 ). 
Any emission with a canonical shell-type spectral index be- 
tween —0.5 and —1 would be far below our detection limit. 
Emission on the scales Slane et al. expected from the X- 
ray brightness distribution would be heavily resolved ex- 
cept on the shortest of our baselines. Modeling the X-ray 

in the Chandra archive. 



Crab-like SNR G21.5— 0.9 at millimeter wavelengths 



3 



emission with 60"-T30" diameter Gaussian, disk and shell 
distributions, we generated uv data for the most compact 
array configuration used. The expected visibility of the 
extended emission ranged from 5 to 50% for these mod- 
els. After tapering our data to a 1' beam, we were unable 
to detect any emission to a (lcr) level of 10 mJy beam -1 
corresponding to a source brightness limit of 20 to 200 
mJy arcmin -1 for the above models. 

4. DISCUSSION 

The virtually identical radial brightness distributions 
above and below the apparent break frequency given by 
Salter et al. (1989a) indicate no change in the spectral 
index with radius. This would appear to rule out any 
change in magnetic field strength, and thus break fre- 
quency with radius and also any variation in diffusion of 
particles throughout the remnant. Alternatively the break 
may actually be above the frequency of our observations. 
We discuss both possibilities below, but first we revisit the 
questions of the age and distance of G21.5— 0.9. 

The most reliable distance estimates for G21.5— 0.9 have 
been minimum distances obtained kinematically from hy- 
drogen absorption measurements. These estimates have 
ranged from 4.4 to 5.5 kpc. The best determination of 
vlsr remains that of Davelaar et al. (1986) (wlsr = 65 km 
s -1 ). We use R = 8.0 kpc (Reid 1993) and the Galac- 
tic rotation model of Brand & Blitz (1993) to obtain a 
minimum distance to G21.5— 0.9 of 4.4 kpc. At that dis- 
tance the major axis of about 100" is 2.1 pc. See Davelaar 
et al. (1986) for a discussion of the maximum distances 
sustainable from a kinematic analysis. 

There is no reliable age estimate for G21.5— 0.9. While it 
is physically smaller than the other Crab-like SNRs stud- 
ied, it is much fainter than the Crab and also has a low 
frequency for the break in the spectrum, implying that the 
pulsar either did not provide as much energy, or has had 
time to decay. G21.5— 0.9 may be the result of a super- 
nova in a very dense part of the Galaxy, or its progenitor 
may have had a relatively weak stellar wind. The radio 
size alone indicates that the remnant could be young com- 
pared to the historical Crab-like SNRs (the Crab Nebula 
and 3C58). But the larger size of the faint outer X-ray 
component indicates a greater age. A reasonable estimate 
for the age could be 1000 years, although this may be sub- 
stantially in error. 

4.1. Spectral break below 100 GHz 

If the break frequency is below our observing frequency 
and particles are injected locally by the neutron star as 
the nebula expands, we would expect to see some spectral 
steepening at the edge of the remnant, where electrons 
with shorter lifetimes would be under- abundant. This 
would show up as a smaller size, as is seen in X-rays but 
not in the radio. Applying the standard minimum energy 
requirement (cquipartition), the average magnetic field for 
G21.5— 0.9 at our adopted minimum distance is 460y(iG 
(Slanc et al. 2000). This is comparable to estimates of the 
field in the Crab Nebula obtained by this and other meth- 
ods (Hester et al. 1996). The conventional spectral break 
due to synchrotron losses in a field B mG occurs at age 
- 400005~ 1 ' f V-°- 5 yr at frequency v GHz. For v h < 100 

3 Available at http://www.mrao.cam.ac.uk/surveys/snrs. 



GHz the estimated age of the break would be > 13, 000 
yr, larger than seems likely for this small-diameter object, 
unless it was formed in a very dense medium. Magnetic 
fields of several mG are required to make the synchrotron 
lifetime comparable to the age of most known Crab-like 
remnants exhibiting only a naked plcrion. This supports 
the view that this spectral break is above our observing 
frequency. 

Another explanation of the observed spectrum, dis- 
cussed by Green & Scheuer (1992), is that in 3C58 and 
G21.5— 0.9 the pulsar's excitation has stopped and so the 
decay of the older electrons has proceeded further and 
caused the steeper high frequency spectrum. If the rel- 
ativistic electrons do diffuse throughout the remnant, this 
could allow the same structure above and below the break 
frequency A related explanation for the low frequency 
of the spectral break and the steeper spectrum above the 
break comes from a suggestion that Crab-like SNRs might 
evolve to be more like G21.5— 0.9 when the pulsar input 
decays; young electrons observed in X-rays experience a 
weaker magnetic field when they are injected into the pul- 
sar wind nebula than did the older ones still seen at radio 
wavelengths (Woltjer et al. 1997). We note that 3C58 and 
G21.5— 0.9 both have lower break frequencies and steeper 
spectral indices above the break than remnants with active 
pulsars such as the Crab Nebula and 0540 — 693 (Green 
2000 3 , Manchester et al. 1993). 

In summary, if the break is below our observing fre- 
quency, and due to synchrotron losses, we would expect the 
spectral index to steepen with distance from the remnant's 
center. Since this is not observed, either the electrons were 
not injected with a simple power law and perhaps have an 
intrinsic spectral flattening towards longer wavelengths, or 
the radio-emitting particles have diffused throughout the 
remnant. 



4.2. Spectral break above 100 GHz 

The above discussion and the X-ray results seem to fa- 
vor the conclusion that the break is indeed above approx- 
imately 100 GHz. Perhaps the slightly lower integrated 
flux densities in the 50-100 GHz range could be accounted 
for by different background levels chosen by different ob- 
servers, particularly if there is any faint radio counterpart 
corresponding to the extended X-ray emission. 

Gallant & Tuffs (1998) have measured the spectrum of 
G21.5— 0.9 at infrared wavelengths and find evidence for 
an X-ray to radio change in spectral index of —1.0, which 
could indicate a break frequency much higher than previ- 
ously thought. Extrapolation of the X-ray spectra (Slane 
et al. 2000; Warwick et al. 2001) to meet the extrapolation 
of a single power law fit to all the radio data or just to the 
points below 30 GHz tends to support this conclusion. We 
note that a recent report of an integrated flux density of 
3.9 Jy for the entire SNR at a frequency of 230 GHz by 
Bandiera et al. (2001) also supports a single power law for 
the radio spectrum at least beyond that frequency. Sub- 
millimeter and further infra-red observations are required 
to settle the question. 

Overall, the observational data better support a higher 
break frequency due to simple synchrotron losses, with 
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spatial effects seen in X-rays but not at radio wavelengths. 
In this case, the particles on the outer edges of the rem- 
nant should have been injected/accelerated first, and are 
the ones showing evidence that the highest energy parti- 
cles have begun to decay. 

4.3. The X-ray morphology 

The radio data do not address the question of why there 
are distinct breaks in the X-ray morphology between the 
central peak and the core and again between the core and 
the plateau. The former characteristic is similar to that 
seen in 3C58 (Helfand et al. 1995) but the latter appears to 
be a unique property of G21.5— 0.9. One possibility is that 
there could have been multiple events or injection epochs. 
Although the entire remnant seems to be pulsar-powered, 
expansion into different ISM may affect the way the ma- 
terial responds to the pulsar stimulation. Alternatively, 
the pulsar may have a nebula or a companion which could 
change its effect on its surroundings. The breaks in the 
X-ray spectrum may also be due to different decay rates 
for the emitting particles in the core and the plateau. The 
recent X-ray data argue against any thermal component 
(Safi-Harb et al. 2001). 



5. CONCLUSION 

Using BIMA observations at 94 GHz we have found 
G21.5— 0.9 morphologically to have very similar structure 
at millimeter and longer radio wavelengths, indicating no 
change in spectral index with radius (although we have not 
been able to reproduce all the small-scale structure seen at 
22 GHz). Particles emitting at radio wavelengths either do 
not show evidence of a variation in lifetime with distance 
from the assumed pulsar, or are well-diffused through- 
out the remnant. Since radial variations in the spectrum 
are seen at X-ray wavelengths, the most likely scenario 
is that the break frequency is above the frequency of the 
present observations. Measurements at higher frequencies 
will most likely support this picture. If they instead con- 
firm a break frequency well below 94 GHz, the break is 
not likely to be due to the limited lifetimes of the higher 
energy electrons, but to some other mechanism. 
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Fig. 1. — 94 GHz mosaic image of G21.5— 0.9 (contours) overlaid on (a) greyscale of 94 GHz data and (b) greyscale of 22 GHz Nobcyama 
data. The contour interval is 7 mjy beam -1 , with the lower contour at 7 mjy beam -1 . The 22 GHz image has been normalized to the BIMA 
intensity scale (see text). The restoring beams for each image are shown. 
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Fig. 2.— Radial profiles of G21. 5-0.9 at 22 GHz (Nobcyama), 94 GHz (BIMA), and in the 0.1-10 kcV band (Chandra HRC). The 94-GHz 
symbol sizes represent the uncertainty due to imaging errors. 



